Quantum percolation in quantum spin Hall antidot systems 
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We study the influences of antidot-induced bound states on transport properties of two- 
dimensional quantum spin Hall insulators. The bound states are found able to induce quantum 
percolation in the originally insulating bulk. At some critical antidot densities, the quantum spin 
Hall phase can be completely destroyed due to the maximum quantum percolation. For systems 
with periodic boundaries, the maximum quantum percolation between the bound states creates 
intermediate extended states in the bulk which is originally gapped and insulating. The antidot in- 
duced bound states plays the same role as the magnetic field in the quantum Hall effect, both makes 
electrons go into cyclotron motions. We also draw an analogy between the quantum percolation 
phenomena in this system and that in the network models of quantum Hall effect. 
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Quantum spin Hall (QSH) insulators have been known 
as two-dimensional (2D) topological insulators, which 
have been studied intensively because of the novel physics 
they hosti - — The bulk of these systems shows insulating 
gaps but the system boundaries processes time-reversal 
symmetry (TRS) protected gapless surface or edge 
modes, characterized by the Z2 topological invariant. 4 9 
Disorders are believed not able to cause backscattering 
of these modes because of their spin-momentum locking. 
In two dimension, these gapless modes are also known 
as helical edge states that are responsible for the QSH 
effect 4 6 . In our previous works, we have shown that 
defects in the form of vacancies in 2D topological insula- 
tors can induce bound states that appear within the bulk 
band ga p 1Q i n . Electrons are able to go into cyclotron mo- 
tions around the vacancies through these bound states, 
which is caused by the non-trivial underlying topology 
of the system. In reality, defects inevitably exist in 2D 
QSH systems. It is however not particularly clear how 
these bound states affect the QSH phase as well as the 
transport properties of the helical edge states. 

In this paper, we study QSH systems in which vacan- 
cies appear as arrays using theoretical models. Such kind 
of man-made systems have been referred to as antidot 
arrays i 12 i 13 In the well-studied 2D electron systems such 
as GaAs/GaAs based and InAs based hetero junctions, 
antidot arrays display intriguing physical properties* 1 ^ 
QSH insulators are ideal platforms for studying anti- 
dot systems since they are true 2D materials with ex- 
otic physical properties. By defining antidot lattices on 
the 2D QSH systems, we show that antidots induce ex- 
tra energy bands which can greatly smear the band gap. 
When exceeding the critical antidot density, the QSH 
phase can be completely destroyed. When antidots are 
randomly distributed, the existence of antidot-induced 
energy bands and band gap smearing is confirmed by 
calculating the density of states of the 2D system. We 
also study the transport properties in these systems un- 
der different boundary conditions. It is found that quan- 
tum percolation occurs taking the bound states as step- 



ping stones. Physically this happens because of the fi- 
nite spatial distribution and overlapping of these bound 
states. In the open boundary condition (OBC), the 
bound states formed transverse percolation channels al- 
low the backscattering between conducting channels at 
opposite edges of the sample with the same spin, which 
then destroys the QSH phase (see Fig. le). In the closed 
boundary condition (CBC), the longitudinal percolation 
forms new conducting channels between the source and 
drain in a two-terminal configuration, as illustrated in 
Fig. Id, which can be viewed as truly extended states. 
We found that transition from QSH phase to conven- 
tional insulating phase accompanies appearance of these 
extended states. 

Following our previous wor k] 10 ! 11 , a modified Dirac 
model is used to describe the 2D QSH insulators 15 , 



H () = vp ■ ex + (mv 2 - Bp 2 ) f3, 



(i) 



where —Bp 2 — —B(p 2 -\-p 2 ) is the quadratic correction to 
the topological mass term mv 2 , pi = —ihdi (i G {x, y}) is 
the momentum operator, v, m and B have the dimension 
of speed, mass and reciprocal mass, respectively. The 
Dirac matrices satisfy the anti-commutation relations 
ai<%j = — otjOLi (i 7^ j), di/3 = — fioii and a 2 = j3 2 = 1. 
A common representation of the Dirac matrices can be 
expressed as, 



OLi cr x 



><7i, P = cr z 



> (To, 



(2) 



where cfi =x ^ y ^ z are the Pauli matrices, <Jo is the 2x2 
unit matrix, and represents the Kronecker product. 
The Hamiltonian (pQ) preserves time-reversal symmetry 
THqT -1 = Ho, where T — —ia x a z K is the time- 
reversal operator and K is the complex conjugate op- 
erator. The topological nature of this Hamiltonian has 
been discussed in detail in our previous works 15 . This 
system is topologically non-trivial (with Z2 = 1 ) while 
mB > 0, and trivial (with Z2 = ) while mB < 0. 4 
After appropriately reordering the basis, the Hamilto- 
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FIG. 1: (Color online) Schematic description of the influences 
of in-gap bound states in 2D TIs. (a) Bound states induced 
by a single antidot. Red (black) colors represent the spin- 
up (-down) component; (b) Antidot lattice with CBC that 
is infinite in x— direction. Additional antidot induced energy 
bands appears in the bulk band gap; (c) Antidot lattice with 
OBC, the antidot bands co-exist with the helical edge states; 
(d) In CBC, at a critical antidot density, the electron perco- 
lates through the overlapping bound states and form longitu- 
dinal conducting channels; (e) In OBC, transverse percolation 
channels allows backscattering between edge states on the op- 
posite edges and finally destroy the QSH phase. 



nian (pp) can be decoupled into two block o 10 ^ 1 , 

h± = (mv 2 - Bp 2 )cr z + v(p x a x ±p y CF y ), (3) 

with h- being the time-reversal counterpart of Be- 
cause the system keeps TRS, we only consider the h+ 
block. When being discretized onto an square-lattice, 
the tight-binding version of h+ reads, 

^ B = E c tt MC 't+E(CM+ Ec -) 

i i 

+ E( C iWT t ya t + H.c.) (4) 

where = (c^,d^) T is the annihilation operator of 
electrons at site i with up spins but belonging to two 
different bands 6 , S x (5 y ) is the unit vector displacement 
(with length a) between two nearest-neighbor sites in the 
x (y) direction, and the on-site energy and hopping ma- 
trices in x (y) directions reads, 



„ / 2 Wh 2 \ 
M= Imv - )°z, t x(y) 



ihv 
~ 1 2a C 



BH 2 



-cr z . 



(5) 



Numerically, the antidot is modeled by leaving the 
antidot regions as vacancies on the tight binding lat- 
tice. The CBC geometry is realized by adopting pe- 
riodic boundary condition in the y direction (see Fig. 
1). Transport features are investigated through the com- 
bination of the non-equilibrium Green's function^ and 
the Landauer-Biittiker 18 formalism. During the calcu- 
lation, the effects of each semi-infinite leads are taken 
into account by introducing the corresponding self-energy 
terms. To clearly illustrate the physics, we adopt such 
model parameters to describe the clean substrate: mv 2 = 
— 10, Hv = 72, and Bh 2 = —28, which is in the topologi- 
cally non-trivial regime, i.e., mB > 0. In the tight bind- 
ing calculations, the lattice constant is set to be a = 1. 

For antidots on a finite size sample, its induced bound 
state can be identified by directly diagonalizing the total 
tight-binding Hamiltonian. The bound states are states 
that have spatially localized wave functions and discrete 
energy spectrum. The antidot-induced bound states in 
2D topological insulators are found both in the bulk band 
and in the band gap with wave function confined around 
it (Fig. 2a) 1Q i n . However, those with energies existing 
in the bulk band gap are of most interest to us since they 
coexist with helical edge states in the QSH insulator. As 
an illustration, we plot out the in-gap bound state wave 
function for a single antidot and a simple antidot array in 
Fig. 2. As is shown, the bound states are indeed localized 
around the antidots. Away from the center their wave 
functions quickly decay, which can be approximated in 
the form ipe~ r ^ The characteristic decay length 
£ in Fig. 2 in theory is about 11a, which agrees well with 
our numerical results. For the antidot arrays, as is seen in 
Fig. 2, the in-gap bound states they induce overlap with 
each other and eventually become connected. 

Consider an antidot array with 100-sites in width. The 
finite-size effect-induced energy gap in the clean system 
of this size is around 1.6 x 10 -3 , 16 , which is negligible 
compared with the bulk gap (2mv 2 = 20). The antidots 
are equally spaced in both x and y directions with the 
spacing being D sites. We calculate the spectrum of the 
antidot lattice with D = 5,4,3,2 under both open and 
closed boundary conditions. Meanwhile the wave func- 
tion of the lowest conduction band at = ir/(D+l) (red 
solid circles in Fig. 3) is taken as an example to demon- 
strate the emergence of the intermediate states formed 
by the bound states in the band gap (Fig. 3 insets). 

As the antidots become dense, the bulk band gap is 
found to close up and then re-opens, as is shown in Fig. 
3. When D = 5 and 4, the edge states survive under 
OBC. Meanwhile the antidot-induced energy bands al- 
ready appear between E = ±10 where the bulk band 
gap was, as is seen in the lowest conduction band wave 
functions in CBC in the inset of Fig.3e. When D grad- 
ually decreases, the band gap shrinks more. At D = 3 
the edge states in OBC disappear and the band gap is 
no longer non-trivial (Fig. 3c). In this case, the lowest 
conduction band becomes a bulk mode, as shown in the 
inset of Fig. 3(c). Since the antidot size cannot be con- 
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FIG. 2: (Color online) Connecting of in-gap bound states as 
the antidots become dense. 2x2 sized antidots are equally 
spaced on a 20 x 72 torus geometry, (a)-(d) show the evolution 
of overlapping for in-gap bound state wave functions with 
antidot number Nd = 1 to 4. 



tinuously changed on the lattice model, this implies that 
there exists a critical point between D = 4 and D = 3 
where the band gap vanishes and the QSH phase is de- 
stroyed. As the antidots become denser, the band gap 
further opens up and no states can be found in the low 
energy regime. 

After consider the physical picture illustrated in Fig. 2, 
it can be inferred that there is a critical point where 
the maximum percolation caused by the antidot-induced 
bound states exists. This can be explicitly demonstrated 
by calculating the two-terminal transmission coefficient 
Ti r with antidots randomly distributed in the sample. 
We choose a 2D sample composed by W x L = 160 x 160 
sites, lxl sized antidots are randomly patterned at a 
density p, which is denned as the number of antidot sites 
divided by WL. For example, in the antidot array case 
of Fig. 3c and Fig.3g, the antidot density is p = 1/16. 
The Fermi level is fixed at E^ w but carefully avoided 
from the tiny gap from finite-size effect. The T\ r versus 
p curves under OBC and CBC are shown in Fig. 4(a) 
and (b), respectively. For each p, the sample averaging is 
taken with 100 random antidot configurations. Indeed, 
when p is small, T\ T remains one/zero under open/closed 
boundary condition. But when p approaches a critical 
value p c ~ 9/160, under OBC T\ r drops dramatically 
while under CBC T\ r acquires a sharp peak with the peak 
value approaching 1, which is a signature of maximum 
percolation in the bulk. When p ^> p c , T\ r vanishes 
under both boundary conditions. The critical point value 
we find p c = 9/160 agrees very well with that we inferred 
from the previous results from the antidot arrays in Fig. 3. 





FIG. 3: (Color online) Energy spectrum of QSH antidot lat- 
tices at different antidot spacing D under OBC and CBC. The 
width is 100-sites and the antidot size is 1 x 1 site, (a)-(d) 
are for OBC, (e)-(h) are for CBC. The solid red circle in each 
figure indicate — 0.1 jj^pn on the lowest conduction band, 
the corresponding wave function for three principal layers are 
shown in the insets. 



And thus the physical picture of bound state induced 
percolation can be confirmed. 

The above picture is further confirmed by calculating 
the system's density of states. To eliminate the contri- 
bution of the edge states in the bulk band gap, CBC in 
both x and y directions are adopted. We choose three 
typical antidot densities: p\ = 5/160, p c = 9/160 and 
P2 = 15/160. For each density, the averaged density 
of states is calculated with 100 random antidot config- 
urations using the kernel polynomial method^. When 
p < Pc the original inverted bulk gap is reduced but still 
survives to protect the QSH phase, as is shown in Fig. 
4(c). For p = p ci a sharp peak appears at E = as is 
shown in Fig. 4(d), which is from the conducting channel 
in the bulk formed by the in-gap bound state. The gap 
seen between the conduction and valence band in Fig. 
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FIG. 4: (Color online) Quantum percolation of electron 
through in-gap bound states induced by randomly distributed 
antidots. (a) and (b) show the two terminal transmission co- 
efficient Ti r v.s. p curves under OBC and CBC with sample 
averaging. The sample is W x L = 160 x 160 in size, (c)-(e) 
are the corresponding density of states at pi = ^§q, Pc — y§q 
and pi — respectively. In these calculations, only the /i+ 
sub-block is considered. When the time-reversal symmetry is 
conserved, the h- block produces the same result. 

4(d) is solely due to the finite size effect, because it can 
be inferred from Fig. 3(g) that in the continuum limit 
both the valence and conduction band spectrum should 
become continuous. When p > p c , the band gap re-opens 
and the system becomes a conventional insulator [Fig. 
4(e)]. 

It is well established in the quantum Hall effect 
that the quantum percolation process is responsible for 
the transition between the quantized Hall conductance 
plateau o 19 i 2Q , which is a direct manifestation of the 
localization-delocalization quantum phase transition. In 
the quantum Hall effect, electrons go into cyclotron mo- 
tions because of the high perpendicular magnetic field in 
which process the disorders help to pin the cyclotron or- 
bitals locally. The quantum percolation can thus be intu- 
itively understood in the sense that the cyclotron orbitals 



connect with each other at some critical external condi- 
tions (e.g.: magnetic field strength, Fermi surface, etc). 
Based on this physical picture, the Chalker-Coddington 
network (CCN) model is developed which gracefully de- 
scribes the localization-delocalization transition in the 
quantum Hall effect in the presence of disorder a 19 i 21 . 
Upon each of such transition, it is also accompanied by 
the change of the topological order of the system, namely 
the Thouless-Kohmoto-Nightingale-Nijs integer 23 , which 
changes by 1 at each of these transitions. 

In the QSH insulator we studied, the maximum perco- 
lation accompanied with destroying of the QSH phase can 
also be seen as such kind of localization-delocalization 
quantum phase transition. But in this case, it is accom- 
panied by the change in the Z2 topological invariant from 
1 to , 24 i 25 Through our previous papers, it is made clear 
that the bound states themselves have helical nature (or 
chiral when considering only h+ or h- sub-block ) 1Q i n ; 
which is similar to the chiral nature of cyclotron orbitals 
in the quantum Hall effect. By analogy it is possible that 
at the critical antidot density, the bound states connect 
and form 2D networks in real space, i.e. the formation 
of percolation. At this condition, truly extended states 
exist and the localization length of electrons should di- 
verge. The percolation provides new conducting channels 
in the bulk, which allows edge states be scattered back- 
ward across the sample through the other edge in OBC 
[as is seen in the drop of T\ r in Fig. 4(a)]. Meanwhile, 
the existence of extended states is proven by the dra- 
matic raise of T\ r shown in Fig. 4(b). It is noted that in 
the QSH insulator, there are always two copies of perco- 
lation channels with opposite chirality, one from /i + and 
the other from h-. 

To conclude, our results indicate that appropriate den- 
sity of antidots in a 2D QSH insulators can dramati- 
cally change the sample's electronic structure and con- 
sequently the transport features. There exists a critical 
antidot density at which the bulk electrons becomes truly 
delocalized because of the quantum percolation channels 
across the sample formed by the antidot induced bound 
states. This localization-delocalization transition is ac- 
companied by the vanishing of the QSH phase. 
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